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Abstract 

Associated production of charged gauge bosons and charged charm¬ 
ed mesons at the LHC is considered in the framework of fej-factorization 
approach . Theoretical predictions are compared with ATLAS data, 
and reasonably good agreement is found. Predictions on the same-sign 
configurations are presented including single parton scatter¬ 
ing and double parton scattering contributions. The latter are shown 
to dominate over the former, thus giving evidence that the proposed 
process can serve as another indicator of double parton interactions. 
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1 Introduction 


Having the LHC put into operation, one got access to a number of ‘rare’ pro¬ 
cesses which would have never been systematically studied at the accelerators 
of previous generations. In this article we draw attention to the associated 
production of weak gauge bosons and charmed mesons. This process is in¬ 
teresting on its own as providing a complex test of perturbative QCD and 
our knowledge of parton distributions. Moreover, we argue that it can serve 
as an indicator of double parton interactions, nowadays widely discussed in 
the literature [11013]. 

This article was greatly stimulated by the recent measurement of the 
WD production cross sections reported by ATLAS Collaboration |1]. In 
that study, the interest was mainly focused on the properties of strange sea 
(see discussion below) and, therefore, in order to suppress other possible con¬ 
tributions (considered in this context as background), the authors have only 
presented the difference between the opposite-sign (OS) and same-sign (SS) 
WD production cross sections, In particular, this excludes the 

Double Parton Scattering (DPS) processes which yield same-sign 
and opposite-sign combinations with equal probability. 

On the contrary, we are more interested in just detecting the DPS events, 
and so, will lay emphasis on the SS states. Our article is organized as follows. 
First, we describe our theoretical approach and check its validity by compar¬ 
ing with ATLAS data on Then we extend our consideration 

to the SS states and make predictions for the production cross sections and 
kinematic observables which could be useful in discriminating the SPS and 
DPS contributions. 


2 Theoretical framework. 

Comparison with ATLAS data. 

At the parton level, the production of opposite-sign states is dom¬ 

inated by the quark-gluon scattering 

g + q ^ W~ -he or g + q ^ W^ -h c (1) 

followed by nonperturbative fragmentation of c-quarks into charmed mesons. 
Here the main contribution comes from strange quarks, while the contribu¬ 
tion from d-quarks is suppressed by Cabibbo angle. 
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To calculate the physical cross sections, we employ the fc^-factorization 
approach HIE]. Here we see certain advantages in the fact that, even with 
the leading-order (LO) matrix elements for hard subprocess, we can include 
a large piece of next-to-leading order (NLO) corrections taking them into 
account in the form of /crdependent parton densities. In this way we auto¬ 
matically incorporate the initial state radiation effects, which play important 
role in the event kinematics. Further on, the formally NLO subprocess 

g + g^W~+c+s or g + g^W~^+c+s (2) 

no longer needs to be added because it is aleady contained in (1). Indeed, 
the quark-gluon coupling in subprocess (2) can also be regarded as part of 
the evolution of sea quark densities q{x) and q{x) in (1). 

On the technical side, our calculations follow standard QCD and elec- 
troweak theory Feynman rules, but the initial gluon spin density matrix is 
taken in the form HI M = k^k^/\kT\^, where k^ is the component of 

the gluon momentum perpendicular to the beam axis. In the collinear limit, 
when /ct’ —)■ 0, this expression converges to the ordinary ege*^ = — | g^^’^, 
while in the case of off-shell gluons it contains an admixture of longitudinal 
polarization. 

In our numerical analysis we used KMR |7] parametrization for uninte¬ 
grated quark and gluon distributions with MSTW [8] collinear densities taken 
as input; we used running strong and electroweak coupling constants normal¬ 
ized to Q;<j(m|)=0.118; a(m|)=l/128; sin^ ©vc = 0.2312; the factorization 
and renormalization scales were chosen as = rri^+p^{W)] 

the c-quark mass was set to mc=1.5 GeV; c-quarks were converted into D^*'^ 
mesons using Peterson fragmentation function [9] with e = 0.06 and normal¬ 
ized to /(c D) = 0.268 and /(c —)■ D*) = 0.229 [TU] . 

Our results obtained with this parameter setting are summarized in Table 
1, where we present the production cross sections integrated over 

the hducial phase space region described in Ref. [1]. We observe reasonable 
agreement with ATLAS data. 

3 Same-sign states 

and double parton interactions 

Now, having our approach validated, we turn to double parton scattering. 
Detecting same-sign conhgurations is certainly preferable here, be- 
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Table 1: Measured and predicted cross sections times the W branching 
ratio (in pb) integrated over the hducial region Pt{1) > 20 GeV, \ri{l)\ < 2.5, 
Priiy) > 25 GeV, Pt{D^*^) > 8 GeV, < 2.2. 



Data 

Theory 

BrW^^‘'aO^-sS(^W+D-) 

17.8 

17.7 

BrW^i^aO^-sS(w-i)+) 

22.4 

19.5 

BrW^i^aO^-sS(w"^ D*-) 

21.2 

15.1 

BrW^i^aO^-sS(w~D*^) 

22.1 

16.8 


cause we are then free from SPS background due to subprocesses (1) or (2). 
There are, however, still many other background sources, both direct and in¬ 
direct. Among the direct ones, we consider the quark-antiquark annihilation 
at 0{a‘^a) 


u + d ^ hh^-l-c-l-c or d + u ^ W +c+c (3) 

and quark-gluon scatering at 0{ala) 

g+u —>■ W^+d+c+c or g+d W~+u+c+c. (4) 

Subprocess (jlj) has one extra as in comparison with ([3]), but it employs gluons 
which are more abundant than antiquarks in the proton, and that is why may 
take over. Among the indirect sources we have gluon-gluon fusion 

g + g^W~+c+b or g + g^W~^+b+c (5) 

followed by the decays b —)■ c-|-X or 6 —>• c+X, and the production of top 
quark pairs 

g + g ^ t + t and q + q ^ t + t (6) 

followed by a long chain of decays, such as f —?• W^+b, c-|-s, b —>■ c-fX 

or 6 —>■ c+c+s (and the charge conjugated modes). Here same-sign 
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configurations may be formed by a W~^ boson coming from t and a c-quark 
coming from b coming from t, or a c-quark coming from b coming from i. 
Similarly, in the case of single top production 

u + d ^ t+ b oi d + u ^ i + b (7) 

the configuration may be formed by a boson coming from t and 

a c-quark coming from b coming from the same t, or a c-quark coming from 
b. Note that the subprocesses ([6]) are purely strong, and so, may have large 
cross sections in spite of large f-quark mass. All other possible processes 
beyond ©-(IZl) are expected to be suppressed by extra powers of coupling 
constants (already the case of ([7])) or by Kobayashi-Maskawa mixing matrix 
(already the case of ([5])). 

A comment is needed on the choice of renormalization scale in ([3]). This 
process factorises into the production of W+g* at fi^=m‘^(W) and the sub¬ 
sequent gluon splitting g*^cc, for which the cc invariant mass seems to be 
a more suitable measure. Note that using different values for these two 
different steps does not violate the overall gauge invariance. So, we calcu¬ 
late the resulting cross section with as{rn^{W))as{mlg), regarding it as the 
pessimistic (the upper) limit for the background. By the same reasoning, we 
adopt a^{rri^{W))as{ml^) for subprocess (jl]). 

For the hducial phase space of Ref. [1], we estimate the above contribu¬ 
tions to W~^D~^ and W~D~ states as 


= 0.41 pb; 

(8) 

{du^Wcc) = 0.29 pb; 

(9) 

Br^^^^a^^^^igu^Wdcc) = 1.0 pb; 

(10) 

{gd^Wucc) = 0.7 pb; 

(11) 

{gg^Wbc) = 0.002 pb; 

(12) 

{gg^Wbc) = 0.002 pb; 

(13) 

Br^^^-a^^^'^igg^ti) = 1.1 pb; 

(14) 

{gg^tf) = 1.1 pb; 

(15) 

{qq^tt) = 0.6 pb; 

(16) 

{qq^tt) = 0.6 pb; 

(17) 

Br^^^'^a^^^^iud^tb) = 0.06 pb; 

(18) 

{du^bB) = 0.04 pb. 

(19) 
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The results (]5 |) - ffTTD were obtained assuming the already mentioned frag¬ 
mentation probability /(c D) = 0.268. The results were ob¬ 

tained under the assumption of 100% branching fraction for t —)■ bW, of 
equal fragmentation probabilities for b ^ and b —)■ B~, and using the in¬ 
clusive branching fractions Br{B^ —)■ D^X) = 37%, Br{B^ —)■ D^X) = 3%, 
Br{B~ —)■ D^X) = 10%, Br{B^ —)■ D^X) = 2.5% listed in the Particle 
Data Book [m. The quark masses were set to mt=175 GeV and m6=4.8 GeV. 
We make no predictions for W^D*^ states for the reason of not knowing the 
relevant B —)■ D*X decay branchings. Variations in /i|. and within a fac¬ 
tor of 2 around the default value make a factor of 1.6 increasing or decreasing 
effect on the estimated production rate. However, these effects mostly cancel 
out in the signal to background ratio. 

Now we proceed to discussing the expected signal from double parton 
interactions. Under the hypothesis of having two independent hard partonic 
subprocesses A and i? in a single pp collision, and under further assumption 
that the longitudinal and transverse components of generalized parton dis¬ 
tributions factorize from each other, the inclusive DPS cross section reads 
(for details see, e.g., the recent review [T] and references therein) 


a 


AB _ 
DPS ~ 


A B 

m Usps'^sps 

2 Ueff 


( 20 ) 


where Ueff is a normalization cross section that encodes all “DPS unknowns” 
into a single parameter which can be experimentally mesured. One can 
identify aes with the inverse of the proton overlap functions squared: 


^efF — 




( 21 ) 


where T(b) = J /(bi)/(bi—b) d^bi is the overlap function that characterizes 
the transverse area occupied by the interacting partons, and /(b) is supposed 
to be a universal function of the impact parameter b for all kinds of partons 
with its normalization hxed as 

j /(bi)/(bi-b) d% d% = j T(b) d% = 1. (22) 

A numerical value of aes — 15 mb has been obtained empirically from hts to 
pp and pp data [121131 [H na [T6] . It will be used in our further analysis, 
although an estimate as low as cXeff — 5 mb is also present in d- 
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The inclusive SPS cross sections a^pg and crfpg for the individual partonic 
subrocesses A and B can be calculated in a usual way using the ordinary par- 
ton distribution functions. The symmetry factor m equals to 1 for identical 
subprocesses and 2 for the differing ones. 

In our present case, the inclusive production cross sections a{D^) and 
a(W^) have been calculated in accordance with Refs. [lU] and [IS], respec¬ 
tively. For the considered hducial phase space our expectations read 


^incl(^^) = ^incl(^ ) = 11-4 hb, (23) 

Vinci(hF+) = 3.5 nb, (24) 

Rr^^'Vinci(hF“) = 2.5 nb, (25) 

where the estimates fl24p - fl25p are supported by direct recent measurement 
[T9] , and so, 

Rr^^^VflP5(TF+R+) = 2.7 pb (26) 

Br^^^^aopsiW-D-) = 1.9 pb. (27) 


These numbers are close to the combined SPS contribution. This means 
that the excess brought by DPS to the visible cross-sections is not 

large enough to unambiguously testify for its presence: the DPS signal is 
large, but the background uncertainties are also large. Moreover, the shapes 
of the DPS and SPS kinematic distributions are rather similar: the decays 
of heavy t-quarks and W bosons make the hnal state distributions broad 
and smooth. Selection cuts on the azimuthal angle difference A0 or rapidity 
difference A?/, so promising in other reactions [20], remain practically useless 
in the present case. This is illustrated in Fig. 1 where we show correlations 
between mesons and muons coming from bosons in same-sign events 
at ATLAS conditions. 

Fortunately, the indirect contributions can be signihcantly reduced (if 
not rejected completely) using a well-known experimental technique based 
on the property that the secondary 6-decay vertex is displaced with respect 
to the primary interaction vertex. We are then left with direct background 
(j8 |) - (fTTD lying well below the DPS level, even with conservative estimate of 
cTefr=15 mb and with ‘pessimistic’ choice of hr as dicussed above. In fact, 
our numbers represent the upper edge of the background uncertainty band. 
A similar relation is seen for D* mesons: 

Br^^^'^aDPs{W+D*+) = 2.3 pb (28) 
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Br^^^^aopsiW-D*-) = 1.6 pb (29) 

= 0.35 pb (30) 

Br'^^^^a^~^*~{du^Wcc) = 0.25 pb (31) 

Br^^^^(j^^^*^{gu^Wdcc) = 0.85 pb (32) 

{gd^Wucc) = 0.60 pb (33) 


We thus come to an important conclusion that the production of same- 
sign W^D*^ states is very indicative as DPS signal. This situation is close 
to the production of same-sign pairs proposed earlier in Ref. [21] . 

However, the events occur at a much lower rate and are then less 

convenient for analysis. 


4 Conclusions 

We have considered the production of a W boson in association with a 
charmed meson in pp collisions at the LHC and made a comparison with 
experimental results. Our theoretical calculations have shown reasonable 
agreement with ATLAS data on We have extended our con¬ 
sideration to the same-sign configurations and found that after re¬ 

jecting the 6-decays the DPS signal clearly dominates over SPS background. 
Thus, we come to an important conclusion that the production of same-sign 
states can serve as a new reliable indicator of double parton scatter¬ 
ing. 
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Figure 1: Kinematic correlations between muons and D-mesons in same-sign 
events at ATLAS conditions: distributions in the azimuthal angle 

difference A0 (upper panel) and rapidity difference A?/ (lower panel). The 
different contributions are represented by: solid curve, double parton scat¬ 
tering; upper dotted curve, gg^ti; lower dotted curve, qq^ti] dashed curve, 
ud^Wcc and du^Wcc; dash-dotted curve, gu^Wdcc and gd^Wucc. 
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